Magnetic interference patterns in superconducting junctions: 
Effects of anharmonic current-phase relations 
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A microscopic theory of the magnetic-field modulation of critical currents is developed for plane 
Josephson junctions with anharmonic current-phase relations. The results obtained allow examin- 
ing temperature-dependent deviations of the modulation from the conventional interference pattern. 
For tunneling through localized states in symmetric short junctions with a pronounced anharmonic 
behavior, the deviations are obtained and shown to depend on distribution of channel transparen- 
cies. For constant transparency the deviations vanish not only near T c , but also at T = 0. If 
Dorokhov bimodal distribution for transparency eigenvalues holds, the averaged deviation increases 
with decreasing temperature and takes its maximum at T — 0. 
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I. INTRODUCTION 

Magnetic interference patterns in superconducting 
junctions originate from quantum coherence of the super- 
conducting state under the applied magnetic field. They 
attract considerable experimental attention and under- 
lie effective studies of various problems of superconduc- 
tivity (see for example^— ). Yet the corresponding mi- 
croscopic theory for short junctions is lacking and the 
present understanding of the results is based mainly on 
the Ginzburg-Landau approach and the tunneling limit. 
A microscopic extension of the results to the low tem- 
perature region T -C T c is required, since the Josephson 
current is to a great extent controlled by discrete An- 
dreev bound states, which are not resolved within the 
Ginzburg-Landau theory. Effects of finite transparencies 
of the transport channels are intrinsically connected to 
the contributions of higher harmonics of the supcrcurrent 
to the modulation, which, therefore, can be adequately 
described only beyond the tunneling approximation. Due 
to the absence of a corresponding microscopic theory, ex- 
perimental data on the interference patterns are analyzed 
in the literature partly phenomenologically with reference 
to usual procedure firmly confirmed for tunnel junctions 
near T c . 

When a transparency of plane junctions gets close 
to unity, there is usually a crossover from the Joseph- 
son current to bulk superconducting flow. Nonetheless, 
there are important plane contacts, where the physics 
of weak links is still valid even in the presence of 
highly transparent transport channels. This is the case 
for long superconductor-normal metal-superconductor 
(SNS) fully transparent junctions of various geometries, 
where interference patterns have been studied in detail 
theoretically at arbitrary temperatures and beyond the 
tunneling approximation^— . In particular, the central 
Fraunhofer peak in clean planar and long SNS junctions 
with fully transparent interfaces has been found to get 
strongly distorted at low temperatures. At T = it 
acquires a triangular for m 11 ! 12 , which correlates with a 
saw-toothed current-phase relation taking place under 



the same conditions in the system o 12 ' 20 . This exam- 
ple demonstrates that pronounced anharmonic current- 
phase relations in superconducting junctions can entail 
significant qualitative modifications in the corresponding 
magnetic interference patterns. 



Another characteristic weak link with a strongly anhar- 
monic current-phase relation is a short clean highly trans- 
parent point contact, which in a fully transparent case re- 
duces to the Kulik-Omclyanchuk clean superconducting 
constriction^!. Similar results also occur for tunneling 
through a single localized state or for plane junctions, 
where resonant electron tunneling takes place via indi- 
vidual localized states homogeneously distributed over an 
insulating interface (see^^r— and references therein). In 
such systems an analytical description of the Josephson 
current is possible at low densities of the transport chan- 
nels with arbitrary transparencies since the pair breaking 
effects are small there. 



In the present paper modulations of the critical current 
are described based on a microscopic theory of Josephson 
junctions generalized to the case of an applied magnetic 
field. An integration of the modulated current over the 
plane of a rectangular junction is carried out explicitly in 
a general form for arbitrary interface transparencies. The 
answer is related to the phase dependent part of the ther- 
modynamic potential in the absence of the modulation, 
taken at the field-dependent phase difference. The theory 
is applied to short junctions with localized states homo- 
geneously distributed over the interface plane. A qualita- 
tive difference is demonstrated between deviations of the 
modulated critical current from the Fraunhofer pattern 
in junctions with constant transparency and Dorokhov 
bimodal distribution of transparency eigenvalues. Junc- 
tions between isotropic s-wave superconductors are con- 
sidered below but the extension of basic results to uncon- 
ventional superconductors is straightforward. 
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II. MODULATION OF THE JOSEPHSON 
CURRENT 

Let superconducting electrodes Si and S r be thick com- 
pared to the magnetic penetration depths A;( r ), while let 
the thickness of the interlayer and the junction width 
be much less than the coherence lengths £u r ) and the 
Joscphson penetration length, respectively. One takes 
the x axis perpendicular to the contact plane and the 
magnetic field applied along the z axis:B(x) = B(x)e z 
(see fig. [1]). It is convenient to take the vector potential 
in the form A(x) = A(x)e y , div A = 0, which coincides 
with the gauge usually taken in describing the Mcissncr 
effect. In contrast to the case of the Meissner effect, in 
Joscphson junctions the vector potential A(x) = A(x)e y 
does not vanish everywhere in the depth of supercon- 
ductors, where the screening supercurrent j y {x) and the 
screened field B(x) do vanish. Indeed, a difference be- 
tween asymptotic values of the vector potential is asso- 
ciated with the magnetic flux $ through the junction: 
A +OQ — A-ca = &/L y . Here L y is a contact width along 
the y axis. 

Nonzero asymptotic values of the vector potential 
can be excluded from microscopic equations by means 
of the corresponding gauge-like transformation. Thus, 
Bogoliubov amplitudes and order parameters can be 
represented as u r ^ = u r( ^ exp [H^4±ooy] , v r W = 
v r & exp [- f c A ±ooV ] , A<"W = A r » exp [^yA ±00 ] . 
Here quantities u r ^\ v r ^ and A r ^ satisfy the equa- 
tions where only the residual parts of the vector po- 
tential A r ( e \x) = A(x) — A± OQ are present in the right 
and the left superconducting regions, respectively. The 
phases of the corresponding order parameters A r W = 
|A r W| exp(ix rW ), A r W = \A r ^\exp(i x r{e) ) are related 
as X r (y) = x r ^ + fr^iooJ/- For nonzero magnetic 
flux one gets A +QO ^ A_oo and the transformation does 
not reduce fully to fixing a gauge since it differs in the 
two regions. Therefore, after excluding constant asymp- 
totic values of the vector potential from microscopic 
equations, the quantities A±oo enter not only gauge- 
dependent phases of order parameters, of Bogoliubov 
amplitudes and of Green's functions, but also a num- 
ber of gauge-invariant physical quantities. In particular, 
as the result of matching corresponding solutions at the 
interface, the phase difference x(y) = X l {y) ~ X r {v) = 
X e ~ X r + i^r {A+oo — A^oo) y enters a secular equation 
and influences the periodic phase-dependent spectrum of 
interface Andreev states. 

After performing the transformation in superconduct- 
ing regions the problem becomes formally more close to 
that of the Meissner effect, since the residual parts of the 
vector potential A r ^ (x) vanish in the depth of the super- 
conductors together with B(x) and j y {x). Similar to the 

problem of the Meissner effect, A(x) in the given gauge 
does not lead to any additional changes of phases of the 
order parameters, even in a strongly nonlinear regime 2 ^. 
For this reason the modulation of the critical current in 
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FIG. 1: Schematic diagram of the junction. 



the microscopic theory is controlled by the spatial depen- 
dence of the phase difference 

x(y) = X + 27r(j,/L,)($/$o), (I) 

where x — X £ ~ X r an d $o = Trhc/\e\ is the supercon- 
ducting flux quantum. 

As the modulation period L y = tt£ b /X = $ o /S(0)A is 
of a macroscopic scale, the quasiclassical theory of super- 
conductivity applies to a microscopic study of the prob- 
lem. Here Ib — (/ic/|eB(0)|) 1 '' 2 is the magnetic length 
and A = d + A; + A r , where d is the interlayer thickness. 
Within the quasiclassical approximation, interface An- 
dreev bound states are associated with coupled incom- 
ing, reflected and transmitted trajectories, which cross 
the interface at one and the same point. In the absence 
of the field, Andreev bound states are degenerate with re- 
spect to the coordinates (yo,Zo) of the reflection points, 
where parallel incoming trajectories with given Fermi ve- 
locity v / cross the junction plane. The total supercurrent 
represents a sum of separate contributions with various 
possible v/. When the external magnetic field is present, 
the quasiclassical boundary conditions, locally applied at 
each crossing point, result in lifting the degeneracy due 
to yo"d e P en clence of the phase difference x(Uo) across the 
interface. The periodic dependence of the quasiparticle 
spectrum on the coordinate j/o of the crossing point, is 
the microscopic origin of the magnetic-field modulation 
of the current. For describing the modulation, one should 
sum (integrate) over yo the contributions to the current 
from respective parallel trajectories for each given vj. 

In the absence of the modulation, the phase-dependent 
part of the thermodynamic potential of the junction can 
be represented as the following sum over Matsubara fre- 
quencies 

oo 

n (T, X ) = -(T/2) J2 toD(ie n , X ). (2) 

n— — oo 

The quantity D(ie„,x) enters the secular equation 
D(e, x) = for eigenenergies of the system and can be 
defined unambiguousl y 26 ! 27 . In the presence of spin de- 
generacy one gets D(e, x) = Q = 2f2 a . 

A variation of the thermodynamic potential with the 
phase difference for a junction under the applied field 
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5f2(T, x, is expressed via the variation Sf2o(T, X ) in 
the absence of the modulation: 



1 
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(3) 

Here a rectangular plane junction is supposed to occupy 
the space (a — L y /2, a + L y /2) along y axis. The param- 
eter a determines a position of the interference pattern 
relative to the junction edges. Since the Joscphson cur- 
rent and thermodynamic potential satisfy the relation 



-2e d 
h dx 



(4) 



the integration of the current over y can be explicitly 
carried out. One obtains from eqs. ^ and (Q| 
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(5) 



Here Xe = X + ^fr - 1S tne en?ec tive phase difference. 
As this follows from eqs. ([2]) and (O, the magnetic-field 
modulation of the Joscphson current at arbitrary temper- 
atures and transparencies is described by the expression 
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(6) 

Eq. © allows calculations of magnetic-field modulations 
of critical currents, provided that the secular function 
D(ie n , x) is known for the junction in the absence of the 
modulation. The secular function can take complex val- 
ues and its property D(—ie n ,x) = E>*{ie n ,x) ensures 
real values of thermodynamic potentials and the current. 
Since eqs. (HJ)- (J3J underly the derivation of Eq. (jB]) and 
have quite general character, eq. ([5]) applies to a vari- 
ety of planar rectangular junctions with any interfaces, 
including those between unconventional superconductors 
and/or with magnetic interlayers. 

In symmetric junctions the Josephson current is car- 
ried solely by subgap states, for which 

N 

6Q (T,x) = 5{-T^ln[2cosh(^( X )/2T)]}. (7) 

i=l 

Here the sum is taken over Andreev state energies 
Ei( X ) > of N transport channels, which can depend 
on trajectory directions and spin indices. According to 
eqs. © and Q, 
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Within its application domain eq. (J5J agrees with eq. (|6|) . 
In particular, eq. (JSJ) reduces to eq. © in the simplest 

case, when D a (ie n ,x) = T\iLi A i [ £ l + E Kx)] and A; 
are independent of X - 

A phase difference Xe,c{T, $), which cor- 
responds to the modulated critical current 
I C (T, $) = \I(T, Xe, c (T, $), $)|, satisfies the equation 

Jo (T, X e,c(T, *) + ^) = Jo (T, Xe , c (T, *) - ^) , (9) 

where Iq(T, x) is the Josephson current in the absence of 
the modulation. In the zero-field limit one obtains from 
eq. (|6]) or eq. (© familiar general relations between the 
Josephson current and the secular function or the spec- 
trum of interface Andreev bound state s 26 ' 27 . It follows 
from eqs. © or (JSJ) that the current always vanishes un- 
der the condition D (ie n , Xe - if) = D (ie n , Xe + if) 

or Ei (xe ~if) = E i (xe + if) ■ Hence, a 27r-pcriodic 
phase-dependent spectrum ensures positions of nodes of 
the modulated Josephson current at $ = n$o, n = 
±1,±2, irrespective of the phase difference. Since 
all even harmonics also vanish at $ = 2n + 1 $o, for such 
values of the magnetic flux the current is formed only 
by contributions from odd harmonics. For small devi- 
ations (5$ of the magnetic flux from rt$0i the current 
and, in particular, its derivative with respect to the 
phase difference always have opposite signs above and 
below each of the nodes. This signifies that the posi- 
tions of minima of thermodynamic potentials as func- 
tions of the phase difference abruptly change by tt at 
<3> = n^Q. Therefore, continuous 0-ir transitions of the 
interference origin take place with varying magnetic flux 
through points $ = n&o (n = ±1,±2,...), where all 
harmonics of the current vanish simultaneously. The 0-tt 
transitions are known to take place, in particular, with 
varying temperature or interface thickness in junctions 
with magnetic interlayers. One can see that such transi- 
tions also take place with varying magnetic field through 
a junction. If the magnetic field, satisfying the relation 
n$o + l)3>o ; is applied, then originally (tt) 

junctions cither evolve to the (tt) state with respect to 
Xe (for n = 0, ±2, ±4 . . . ), or turn into respective tt (0) 
junctions (for n = ±I,±3,±5, ...). This concerns, in 
particular, the standard situation, when the Fraunhofcr 
pattern describes the modulation. 



III. TUNNELING THROUGH LOCALIZED 
STATES 

Consider further nonmagnetic short junctions between 
identical s-wave superconductors, where tunneling via lo- 
calized states with a large broadening occurs. The in- 
fluence of the screening current and the magnetic or- 
bital effects on the Josephson current is usually negli- 
gibly small in such systems, so that the residual vec- 
tor potential can be disregarded. Then the spectrum of 
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FIG. 2: The Josephson current $)/J c ($ = 0) as a func- 
tion of the phase difference in fully transparent (T> = 1) sym- 
metric junctions at T = 0. The current is normalized to the 
critical current at zero-field and taken for various values of the 
magnetic flux through the junction: 1. $ = 0, 2. <f> = 0.25$o, 
3. $ = 0.5$ , 4. $ = 0.75$ , 5. $ = 1.5$ . 

spin-degenerate Andreev states takes the form Ei±(x) = 



±|A|y 1 — T>i sin 2 (x/2), which formally coincides with 
the spectrum of superconductor - insulator - supercon- 
ductor point contacts 2 ^. The transparency T>i is de- 
scribed here by the Breit-Wigner resonance function, 
taken at the energy of the i-th localized state d 24 ' 28 . 
The coefficient T>i can take any value between and 1, 
depending on the energy of the state and its position 
x i Q across the interlayer. Near T c the order parameter 
is small and, expanding all functions in eq. ([5]) in powers 
of E + /T Cl one can keep there only the main quadratic 
term. This leads to the relation I C (T, $) = I c f(T, <£), 
where I c f(T, <E>) describes the Fraunhofer pattern for the 
critical current 



I cF (T,*)=I c {T,0) 



sin 



7T<I> 
$0 



7T$ 
$0 



(10) 



In the particular case I c (T,0)\t->t c = ^E.^>i 
where the sum is taken over possible different vj. 

At low temperatures arguments of hyperbolic func- 
tions in eq. ([U are large. Using the respective asymp- 
totic expressions one obtains within a simplified model 

of constant V: cos Xe,c(0, $) = cos Xc(0, 0) cos (jg^j ■ 

Here the zero-field phase difference is cos \ c (0,0) = 
— (1 — \J\ — T>) 2 /T>. This solution results in the zero- 
temperature critical current, which exactly reduces to the 
Fraunhofer pattern eq. (|10|) for any field value. The zero- 
field critical current at T = 0, which enters eq. (fT0|) as 
a factor and depends on T>, is found to take the form 
J c (0,0) = (|eA|/ft)(l-vT= V). 

Figure [2] shows current-phase dependences for vari- 
ous values of the magnetic flux through symmetric fully 
transparent junctions in question with a few identical 
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FIG. 3: Relative deviations 5I C (T) as functions of tempera- 
ture in symmetric junctions, taken for $ = O^^o and various 
transparencies: l.V = 1, 2.V = 0.95, 3. V = 0.9, 4. V = 0.8, 
5.23= 0.5. 

transport channels at zero temperature. The depen- 
dences involve significant contributions from a large num- 
ber of harmonics. Surprisingly, the conventional in- 
terference pattern for the critical current in symmetric 
junctions takes place in this case. Such behavior con- 
trasts to what is known for long superconductor-normal 
metal-superconductor junctions. Based on eq. ([5]). one 

can calculate relative deviations 5I C (T, $) = [I C (T, <&) — 



I cF (T,$))/\I c (T, $)| of the critical current from the 

Fraunhofer values given by eq. (|10[) . The quantity 
SI C (T, <E>) vanishes identically only in the tunneling ap- 
proximation and/or near T c . Fig. [3] displays the deviation 
SI C (T, $) as a function of temperature, for $ = 0.5<i>o 
and various transparency coefficients. At intermediate 
temperatures eq. (|10[) does not apply exactly, but the 
nonmonotonic temperature dependent deviations due to 
higher harmonics are less than few percent and vanish at 
T = 0. 

In asymmetric junctions zero-temperature deviations 
SI C do not vanish, as this follows from eq. ([5]). They are 
shown in fig.[3]as functions of the parameter 7 = |A;/A r |, 
which characterizes the junction asymmetry. Since 5I C 
does not vary with interchanging left and right order pa- 
rameters, one takes 7 > 1. As can be seen, SI C at T = 
is positive and not large, reaching about ten percent at 
7 = 14 and not exceeding eleven percent even at 7 — 30. 



IV. MULTICHANNEL EFFECTS 

In Josephson junctions with one transport channel or 
a few channels with identical transparencies, the phase 
difference Xe,c, which corresponds to the critical current, 
depends on the transparency value due to anharmonic 
current behavior. In a multichannel junction the critical 
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FIG. 4: Zero-temperature deviations 8I C (T = 0) in asymmet- 
ric junctions as functions of 7 = |Ai/A r |, taken for the same 
set of parameters $ and T> as in fig. [jj] 



current arises from an interplay between channels with 
different transparencies and the resulting value of Xe,c 
can differ from those for the respective one-channel junc- 
tions. For this reason the quantity Xe.c in strongly anhar- 
monic junctions can noticeably influence the dependence 
of the critical current on the transparency distribution 
over transport channels. 

A symmetric junction containing only two transport 
channels with different transparencies X>i, T>i represents 
the simplest example for establishing multichannel ef- 
fects. The corresponding critical current can be found 
based on eq. (J5J . Its relative deviations from eq. (fit))) arc 
shown in fig. [5] as functions of temperature, for T>\ = 1 
and various values of T>2 ■ Characteristic qualitative fea- 
tures of the multichannel case follow from a compari- 
son of figs. [5] and [31 Figure [3J and the first curve in 
fig. [5] show that at intermediate temperatures eq. (|T0|) 
slightly overestimates the modulated critical current in 
symmetric junctions with identical channels (constant 
transparency), but the deviations vanish at zero temper- 
ature. Curves 2-5 in fig. [5] demonstrate that an interplay 
between different channels can change such situation. As 
a result, the critical current exceeds the Fraunhofer value 
at low temperatures and, if T>2 is not too small, maximal 
deviations occur at zero temperature. A nonmonotonic 
behavior of deviations at low temperatures taking place 
with the parameter T>i arises from two competing rea- 
sons. If T>2 is not too small and decreases, then the vari- 
ations of Xe,c from its value for constant transparency 
T>i = 1 increases and determines increasing deviations 
of the critical current from the Fraunhofer one. On the 
other hand, for a channel with sufficiently small and de- 
creasing T>2, its contribution to the total critical current 
diminishes and the deviations are more and more domi- 
nated by one channel with T>\ ~ 1, i. e. they decrease at 
zero temperature up to zero at T>2 = 0. 

The curve 6 in fig. [5] describes the critical-current devi- 
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FIG. 5: Relative deviations 5I C as functions of temperature, 
taken at <E> = 0.5$o in symmetric junctions. Solid curves 
correspond to junctions with two channels, where T>\ = 1 
and T>2 takes the following values: \.T>% — 1, 2.T>2 — 0.9, 
3.V 2 = 0.5, 4.0 2 = 0.3, 5.V 2 = 0.1. The dashed curve 
describes deviations in a junction containing one channel with 
T>\ — 1 and ten identical channels with T>2 =0.1. 



ations in a junction containing ten identical channels with 
T>2 = 0.1 and one fully transparent channel (D\ = 1). 
Contributions to the Josephson current from ten chan- 
nels with T>2 = 0.1 and one channel with T>i = 1 are 
both significant and jointly lead to more noticeable devi- 
ations of the critical current. 

For tunneling through broaden localized states the av- 
eraging of Iq (T, x) over bimodal Dorokhov distribution 
of transparency eigenvalues leads to the current through 
dirty constriction o 24 ' 29 . The corresponding thermody- 
namic potential takes the form 

fio(T,x) = 

= (2irh 2 T/e 2 R N )J2 arcsin 2 f\ A\ sin |/ ^ e\ + | A| 2 ) 

£n>0 



(11) 



and the modulated current I(T, x, ■£) is defined by 
cq. ([5]). In this case the critical current exceeds the 
Fraunhofer value, the deviation SI C (T,^>) increases with 
decreasing temperature and takes its maximum at T = 0, 
as is seen in fig. [51 



V. CONCLUDING REMARKS 

Experimental results for numerous short junctions are 
known to show, as a rule, modulations of the standard 
type, if a spatial distribution of the supercurrent den- 
sity is not substantially inhomogeneousi Prominent 
exceptions include combined 0-7T junctions, vicinities of 
0-7r transitions and special intcrface-to-crystal orienta- 
tions of high-tcmpcraturc or other superconductors with 
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FIG. 6: Relative deviations 5I C (T,&) averaged over the 
Dorokhov distribution of channel transparencies in symmetric 
junctions: 1. T = 0, 2. T = 0.1T C , 3. T = 0.2T C , 4. T = 0.3T C , 
5.T = 0.5T C . 

anisotropic pairings ^ 10 ' 30 . The present calculations al- 
low an extension, in particular, to short junctions with 
inter layers possessing a collinear magnetic order and/or 
between unconventional superconductors. The devel- 
oped approach can be also generalized to take account of 
the current-induced magnetic field resulting in Joseph- 
son vortices in wide junctions. These problems will be 



studied further and published elsewhere. 

In conclusion, a microscopic theory of the magnetic- 
field modulation of the critical current in Josephson junc- 
tions has been developed in the present paper. As a gen- 
eralization of basic microscopic results in the absence of 
the magnetic field, the modulated Josephson current is 
explicitly expressed via a secular function or, for symmet- 
ric junctions, via a magnetic field dependent spectrum of 
Andrccv interface states. Temperature-dependent devi- 
ations of the modulated critical current from the Fraun- 
hofer pattern have been found for short junctions with 
tunneling through localized electronic states. The devi- 
ations depend on transparency distribution over trans- 
port channels. For constant transparency the deviations 
vanish not only near T c , but also at T = 0. Such be- 
havior qualitatively differs from what is known for long 
superconductor-normal metal-superconductor junctions. 
Zero-temperature deviations are found to take place in 
junctions between different superconductors and in sym- 
metric junctions containing channels with different trans- 
parencies. If Dorokhov bimodal distribution of trans- 
parency eigenvalues holds, the averaged deviation in- 
creases with decreasing temperature and takes its max- 
imum at T = 0. It is shown that in a number of junc- 
tions with a pronounced anharmonic current behavior, 
the Fraunhofer pattern is only slightly distorted. 

The support of RFBR grant 08-02-00842 is acknowl- 
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